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We use numerical simulations and semi-analytical methods to investigate the stability and the 
interactions of nontopological stationary qball solutions. In the context of a simple model we map the 
parameter sectors of stability for a single qball and verify the result using numerical simulations of 
time evolution. The system of two interacting qballs is also studied in one and two space dimensions. 
We find that the system generically performs breather type oscillations with frequency equal to the 
difference of the internal qball frequencies. This result is shown to be consistent with the form of 
the qball interaction potential. Finally we perform simulations of qball scattering and show that 
the right angle scattering effect observed in topological soliton scattering in two dimensions, persists 
also in the case of qballs where no topologically conserved quantities are present. For relativistic 
collision velocities the qball charge is split into a forward and a right angle scattering component. 
As the collision velocity increases, the forward component gets amplified at the expense of the right 
angle component. 



I. INTRODUCTION 

It is well known that realistic supersymmetric theories 
are associated with a number of scalar fields with various 
charges. These theories and in particular the MSSM al- 
low for baryonic and leptonic nontopological solitons M 
known as qballs pi- They are composed of squarks, slep- 
tons and Higgs scalars p.H and they originate through 
the breakup of an Affleck-Dine condensate carrying a net 
baryon and/or lepton number [pi. Such objects could 
have interesting cosmological consequences. For example 
they could be responsible for both the net baryon number 
of the universe and its dark matter component W . 

The basic properties of qballs have been studied ex- 
tensively in the literature using mainly analytical and 
semi-analytical methods W|. These methods are partic- 
ularly useful in understanding basic properties of qballs 
like existence, small vibrations, stability in certain pa- 
rameter limits (thick ||] or thin wall [gj approximation) 
but they can not be very illuminating in understanding 
more complicated issues like scattering, interactions or 
stability for arbitrary parameters. The goal of this paper 
is to use numerical simulations of qball evolution in one 
and two spatial dimensions, along with semi-analytical 
methods in order to study the stability under small fluc- 
tuations, the interactions and the scattering of qballs. 

The structure of the paper is the following: In the 
next section we introduce a simple model and show the 
existence of stable qball solutions in its context. By con- 
sidering small fluctuations superposed on these solutions 
for various parameters we construct a map showing the 
stability sectors in parameter space. These sectors are 
then verified by performing numerical simulations of time 



evolution for the solutions considered. A virial theorem 
is also derived and its validity is demonstrated numeri- 
cally. In section III we consider a system of two inter- 
acting qballs and show analytically that the interaction 
potential is time dependent and periodic. This type of in- 
teraction is verified by performing numerical simulations 
of qball evolution which shows that the system performs 
breather-type oscillations with the analytically predicted 
frequency. In section IV we use boosted qball configu- 
rations to perform scattering numerical experiments in 
2-1-1 dimensions showing that for low collision velocities 
qballs tend to scatter at right angles like their topological 
counterparts. For relativistic collision velocities we find 
that the qball charge splits into a forward and a right an- 
gle component. The forward component gets amplified 
as the velocity increases. Finally in section V we con- 
clude, summarize and briefly discuss future extensions of 
this work. 



II. STABILITY - VIRIAL THEOREM 

Consider a complex scalar field <f> in l-f-l dimensions 
whose dynamics is determined by the Lagrangian 



/: = ia„$*9''$-t/($). 



where 

Using a rescaling 
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'The animated simulations may be found on-line at the address http://leandros.chem.demokritos.gr/qballs/index.htm] 
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the Lagrangian (hi) gets simplified as follows 
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where 



B = 



5m ^ 



This leads to the field equation 

$-$" + $- |$|$ 4- B|$P$ = 0. 
Using the usual qball ansatz 

with the conserved Noether charge 
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we obtain the field equation for g(x) 



a" + (^2 _ 1)^ + ^2 
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The requirement of finite energy and the asymptotics ob- 
tained from equation (hfl) imply the boundary conditions 



a'(0) = 0, 
cr(oo) = 0. 



(11) 
(12) 



The combination of equations (|l^ and (11-12) are 
identical to the dynamical equations describing the mo- 
tion of a virtual particle that starts at rest at 'time' x = Q 
from some position ctq and moves to cr = at 'time' 
X — > cx) under the influence of the effective potential 

V^ff ^\{u:^ -im^ + \\^^ -\b\^^ (13) 



By considering the form of the effective potential ( |l3| ) 
it becomes clear that in order to have a solution with 
the boundary conditions (11-12) the following conditions 
must be satisfied 

• Symmetry must be broken by the effective poten- 



tial mm (to 



'eff 



1 <0) 



• There must be at least one intersection point of 
Ve//(cr) with the V;//(|<i>|) = axis. 
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FIG. 1. The field magnitude profile j<3>(a;,t)| of a qball 
solution for B = 4/9, uj'^ = 0.51. 

The first condition implies uj^ < 1 while the second 
condition implies that the equation Ve// = has at least 
one real solution different from 0. It is easy to show that 
the constraints imposed on lo by these two conditions can 
be written as 



9B 



<UJ^ <1. 



The energy density of the system (pT) is 



£ = |$P 






which with the qball ansatz (||) becomes 



S^a'- 



(1 



-c.V-f-^ 



B 



(14) 



(15) 



(16) 



Our main goal is to study numerically the dynamics 
and the stability of single and multi-qball configurations 
in the context of the simple model (|l[) . The first step in 
that direction is to solve equation (O) numerically (in a 
parameter region where qball solutions exist), to find the 
single qball profile and then evolve numerically the so- 
lution according to (0) checking total energy and charge 
conservation. 

Fig. |l| shows the numerically obtained qball profile 
for B = 4/9 and w^ — 0.51 using a fourth order Runge- 
Kutta scheme. The evolution of this configuration in time 
when used as initial condition in equation (R) with peri- 
odic boundary conditions has been performed and found 
to correspond to a stable configuration. The total energy 
and charge Q were conserved to within about 3% during 
the evolution. 

Static scalar field configurations in space dimensions 
higher than two are unstable towards collapse because 
both the gradient and the potential terms of the energy 



may be shown to decrease with a rescahng of the config- 
uration corresponding to coUapse M . On the other hand 
the stabihty of the stationary (time dependent) qball con- 
figuration towards coUapse or expansion is a result of the 
different behavior of kinetic (time dependent) and poten- 
tial terms towards spatial rescaling of the field configura- 
tion. This fact may be seen even in 1+1 dimensions by 
expressing the energy as a sum of kinetic, gradient and 
potential terms 
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After a rescaling of the spatial coordinate x 
expression of the total energy becomes 



1 



(18) 
(19) 
(20) 
ax the where 



a 



(21) 



For stability we demand that the energy is minimized for 
a = I which implies that 



dEc, 
da 



O^h+h 



(22) 



while the second derivative is ^^|q=i = 2 {I1+I2) > 0, 
implying energy minimum and therefore stability. Thus 
the stability of the qball configuration towards coordi- 
nate rescaling implies the validity of the virial theorem 
of equation (22) connecting the potential energy with the 



gradient energy and the kinetic energy (for a similar virial 
theorem see Ref. P]). We have confirmed numerically 
the validity of the virial theorem ( p^ ) and the virial ratio 
was found to be zero to within less than 1% for 



I1+I2 



I1+I2+I3 

various values of uj where qball solutions exist and for 

various B in the range 2/9 < B < 1. 

In order to study the stability of the qball configura- 
tion under small fluctuations that conserve the charge Q 
of equation (0) we must consider variations of the energy 
obtained from the density (|l^) and expressed in terms of 

Q. 



E = dx [a'' 
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(23) 



Since the qball is a stationary solution of the field equa- 
tions, the first variation of (g^) with respect to a vanishes 
identically. The second variation may be written as 



FIG. 2. Parameter sectors of existence and stability for 
qball solutions. 
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and (To is the unperturbed qball solution. For stability we 
require 6^E > for all fiuctuations that conserve charge. 
This is equivalent to demanding that the Hermitian op- 
erator O has no negative eigenvalues. In order to find the 
parameter region corresponding to stability (no negative 
eigenvalues) we have solved the differential equation 



Oda = 



(26) 



with initial conditions 5a (0) — 1 and 5a' (0) = using 
a fourth order Runge-Kutta scheme. By varying the pa- 
rameters B, UJ we have identified the line in parameter 
space where the eigenvalue problem (|26| ) has a ground 
state with zero eigenvalue. This line clearly separates 
the parameter sector of stability from the corresponding 
sector where the equation ( |2q ) has negative eigenvalues 
and the qball is unstable. 

Fig. shows part of the parameter space divided in 
three sectors: the sector where stable qball solutions ex- 
ist, the sector where qball solutions exist but they are 
unstable and the sector where no qball solutions exist 
because the condition ( p^ ) is violated. The scale of the 
plot was chosen in order to magnify the sector of insta- 
bility which would otherwise be too small to be seen. 

We have tested and verified the validity of these sec- 
tors by simulating numerically the evolution of an iso- 
lated qball using (H) in the instability and in the stability 
sectors. 
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FIG. 3. Evolution of a qball solution in the instabil- 
ity sector. The evolution of the field magnitude $(a;,t) is 
plotted. The time evolution was performed in the range 
-7.5 < a; < 7.5. 

Fig. ^ shows an example of evolution of a qball in 
the instability sector. On the other hand we have also 
verified stability by performing the time evolution of 
the field amplitude for a qball in the stability sector 
{u}^ = 0.1 B = 2/9 ~ 0.222). The qball was seen 
to remain stable and evolved undistorted in time even 
though the evolution lasted for ten internal frequency 
periods {tmax = 10(— )). Fig. g] shows the same nu- 
merical experiment for parameters in the instability sec- 
tor (uj^ = 0.05 S = 2/9 ~ 0.222) where equation 
( P6| ) has negative eigenvalues. The evolution lasted only 
three internal frequency periods {tmax = 3( — )) but the 
qball configuration got rapidly distorted and decayed into 
plane waves thus verifying the semi-analytical prediction 
for instability. 



III. INTERACTIONS - SCATTERING 



We now consider multi-qball configurations in order to 
study their interactions and their scattering properties 
using mainly numerical simulations of evolution. For a 
two qball initial configuration we use the ansatz 



$(x,t) 



a+e 



iu^t 



(7-6 



(27) 



where a+ = ai{x + xq) and a- = cr2{x — xo) with ai 
(cr2) the field magnitude for a single qball with lu = cj+ 

(lu = CJ_). 

The interaction potential between the two qballs can 
be obtained by subtracting the energy densities e+(x,t), 
e-{x,t) of each non-interacting qball from the energy 
density of the total interacting configuration. After a 
straightforward calculation we obtain the interaction en- 
ergy as 



FIG. 4. The interaction energy of a pair of qballs calcu- 



lated using the expression of Si, 
showing a behavior of the form - 



,t with B = 4/9, 
- cos(2a;i). 
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where F and G are time independent functions of 0(1) 
and 0{a+(7-) respectively. It is easy to see from the field 
equation dlO) that the fields a± decay exponentially at 
infinity and therefore the term proportional to F domi- 
nates in the interaction energy. We therefore expect that 
the two qball system will perform oscillations with char- 
acteristic angular frequency uj^ — lu- while slowly drifting 
due to the effect of the small constant interaction term G. 
Indeed this is what we see in the numerical simulations 
of the time evolution of the system. 

We have considered the ansatz (p7| ) with uj+ — LOi„t 
and U- = —LOint which implies ai{x) = (T2(a;) since the 
internal angular frequencies of field rotation are opposite. 
The interaction energy of this system obtained numeri- 
cally for Wint = 0.9 as described above by using the nu- 
merically calculated <j{x), is time dependent and is shown 
in Fig. |. 

This interaction energy however assumes that the form 
of the initial ansatz is retained during the time evolution 
of the system and should therefore be subject to test by 
numerical simulation. As expected from the expression 
of Eint the time dependence of the ansatz interaction is 
proportional to cos(2wi„(i). 

In order to test if the predicted time dependence of 
the interaction is realized in a realistic system we have 
performed a numerical simulation of the evolution of the 
ansatz (|^ by solving equation (|^) using a leapfrog al- 
gorithm, periodic boundary conditions and initial condi- 
tions based on (27). The parameter values were B — 4/9, 



— 0.9 and the initial qball distance was 2a;o = 10 in 
a lattice of 2d = 30 (these arc dimensionless due to the 



rescaling of the Lagrangian). The evolution of the field 
magnitude |<i>(x,t)| is shown in Fig. where the oscilla- 
tions of the system can be seen clearly. The period and 
the angular frequency of these oscillations can be found 
by plotting the location of the field maxima as a function 
of time. This is shown in Fig. O where we plot the loca- 
tion Xmax{t) of the maximum of |$(a;,i)| for a; > and 
< t < 15. Clearly the location of the qball at a; > 
described by Xmax (t) performs oscillations due to the in- 
teractions with the qball at x < 0. The period of these 
oscillations can easily be seen from Fig. W to be approx- 
imately T ~ 3.5 which is consistent with the anticipated 
result T = ^ ~ 3.2 based on the interaction potential 
Eint i-B- the spatial frequency of qball oscillation Wspace 
is double the internal frequency uji„t of field rotation. 

We have verified this consistency between the analyt- 
ically predicted period of spatial oscillation and the one 
seen in the simulations for several values of internal ro- 
tation frequency uJmt- 

The result is shown in Fig. M where we plot the ob- 
served angular frequency cospace of spatial qball oscilla- 
tion vs double the corresponding frequency tOmt of inter- 
nal field rotation. The data points can be fitted well by 
a straight line of slope unity as anticipated by the above 
described analytical considerations. During all the sim- 
ulations the energy and total charge of the system were 
conserved to within about 3%. The amplitude of the 
spatial oscillations was not found to be constant but the 
system slowly drifted to larger qball separations with a 
rate dependent on the parameter values. This behavior is 
consistent with the form of the interaction potential Sint 
which includes a subdominant time independent term. 
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FIG. 5. The evolution of the field magnitude correspond- 
ing to a pair of qballs calculated using the ansatz (GTI) with 
B — 4/9, o;f„t — 0.9. As time increases the line becomes 
thicker and the dashes larger. The time range for a com- 
plete period of spatial oscillations is < i < 3 which implies 
bJsr, — 2.1 ~ 2uJint to Within 10%. 
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FIG. 6. Spatial oscillations of the maximum of the field 
magnitude due to qball interactions. 
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FIG. 7. The frequency of spatial oscillations of qball posi- 
tion vs the internal frequency of field rotation. 

In order to study the scattering of qballs we need to 
consider boosted qball configurations i.e. qballs moving 
with a velocity v. Starting from a qball field configuration 
^v=o{x,t) = i7(x)e'"* which does not move in space we 
can construct a configuration <l>„=„„(x, i) that describes 
a qball moving with velocity vq by performing a Lorcnz 
transform to the spacetime variables. Let 7 — . ^ 

(we use units where the velocity of light c is unity) be 
the Lorenz factor and let 



x' = j{x — VQt) 

t' = -f{t - Vox) 



(28) 
(29) 



be the boosted spacetime variables. The field configura- 
tion ^{x',t') expressed in terms oi x, t describes a qball 
moving with velocity vq i.e. ^y^yg{x,t) — $(x',i'). The 
initial ansatz for a two qball system prepared for a colli- 
sion process may be written as 



$(x,t) = $i,=_t,o(x- Xo,t)^-^ +^v=vo{x + xo,t)^^. (30) 

An evolved system of this type with vq = 0.2 is shown 
in Fig. 1^ where we have used uji = ijJ2 = VO.Sl and 
B — 4/9. As it can be seen in Fig. ^ the qball collision 
results in the formation of a long lived central qball which 
subsequently decays to qballs similar to the original ones. 
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FIG. 8. Collision of qballs with identical charges and 
«oi = — ^02 = 0.2. 

A more interesting evolution occurs in two space di- 
mensional systems which will be described in the next 
section. There it will be seen that the effect of right angle 
soliton scattering which has been observed in topological 
solitons persists also in the case of the nontopological 
qballs in a generalized form. 
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FIG. 9. The numerically calculated profile of the qball so- 
lutions for B = 4/9 and various frequency values. In the 
vertical axis we plot the charge density q = u a^. Qballs exist 
for 0.5 < w^ < 1 or 0.71 ^ uj < 1. 



The boundary conditions which have to be met are given 
in Eqs. (11) and (12). Eq. ( p3| ) is expected to have lo- 
calized solutions by arguments which are presented in M 
and are similar to those used for the ID model of Sec. IL 

The search for the solutions is here not straightfor- 
ward as for the ID model (Figs. 1,2). In the present 
two dimensional case, solutions are found by a numerical 
shooting method. Fig. shows the profile of the calcu- 
lated qballs for B — 4/9 and for various values of the 
frequency u. We represent the qball profile through the 
charge density g = wcr^ (cf. Eq. (||)). 

The virial theorem mentioned in Sec. II has, for our 
two dimensional theory, the form 



h^h, 



(34) 



IV. DYNAMICS IN TWO SPACE DIMENSIONS 

The model (P is extended to 2-1-1 dimensions by let- 
ting the index /z take the values /z = 0, 1, 2. The details 
of the formalism are similar to those given in Sec. II. 
That is, the rescaled form of the Lagrangian is given by 
Eq. (0) while the field equation reads 



$ - A$ + $ - |$|$ + S|$p$ = 0, 



(31) 



where A = d^ + dy. We look for axially symmetric solu- 
tions of Eq. (pW and make the qball ansatz 



$(p,t) = a(p)e 



:LiiiJt 



where p — ^x"^ + y^ . The profile cr(p) satisfies 



— + (u?- l)cr ArCp- - Bo^ ^ 0. 



(32) 



(33) 



where the symbols are defined as the two dimensional 
analogues of (18) and (20). The virial relation is used 
to check the precision of our numerical calculations. In- 
deed, the solutions that we find by the shooting method 
satisfy the above virial relation to very good precision, 
better than 1%. 

Finally, we note that our shooting method is unable to 
find the qball profile for the whole frequency range (H). 
Numerical errors render the method inapplicable near the 
lower w-bound. For instance, for B = 2/9 we are able to 
find the solutions with 0.035 < cj^ < 1 and for B = 4/9 
we find those in the frequency range 0.52 < uj^ < 1. 
The difficulties with the numerics should be anticipated. 
Their origin can be traced to the theoretical arguments 
of Ig] with respect to the existence of qballs in dimen- 
sions higher than one. It is actually the friction term in 
Eq. (p3[) which is responsible for the numerical difficul- 
ties. 



We further note that a steadily moving qball can be 
found by applying a Lorenz boost to a static one as dis- 
cussed in the previous section. The calculated qballs will 
be used in the following in numerical simulations on a 
two-dimensional lattice. We shall first study their stabil- 
ity. Then we shall study interactions between two qballs, 
in particular we shall perform simulations of scattering. 

The numerical mesh we use is typically 300 x 300 and 
the lattice spacing is 0.3. Such a numerical mesh is ap- 
propriate to accommodate the structures given in Fig. ^ 
and provide good resolution. The time evolution is per- 
formed by a fourth order Runge-Kutta method. 

First, we put a single qball at the center of our numer- 
ical mesh and simulate its time evolution in order to test 
its stability. We use 2/9 < B < 1 and test all the values 
of Lu for which we are able to find the qball profile by 
our shooting method. We find that qballs are stable ev- 
erywhere in the B — to plane and they travel undistorted 
with the given constant velocity. 

In the next set of simulations we discuss the problem of 
interaction between qballs. Specifically, we focus on sim- 
ulations of scattering between two qballs. Such simula- 
tions have been performed and have proved to be fruitful 
for topological solitons. We shall perform all our subse- 
quent numerical simulations using the values B = 4/9, 
uj = ±0.75. 

We consider a 2D generalization of the ansatz (30) 
which represents two qballs set in a head-on collision 
course. We use wi = W2 = 0.75 and a small velocity, 
namely vq — 0.2. Qballs are represented through their 
charge density: 



Zt 



(35) 



In Fig. 



10 



contour plots are given for the charge den- 
sity at three characteristic snapshots of this numerical 
simulation. The first entry of the figure gives the initial 
ansatz. The two qballs are initially at a distance of 20 
units apart so that there is no overlap and interaction 
between them. Subsequently they collide at the origin 
(middle entry) and scatter at right angles (lower entry). 

The scattering scenario produced by the computer sim- 
ulation is quite interesting. Indeed, this dynamical be- 
havior has been found to be a robust feature in a variety 
of models in two space dimensions which have topological 
soliton solutions WM- However, in our case there is no 
topological invariant associated to the qballs, so we need 
to discuss it further. 

The crucial observation is that the underlying dynam- 
ics which induces this type of scattering can be attributed 
solely to the Hamiltonian structure of the model. It is 
related to the conservation laws of our model and in par- 
ticular to the linear momentum conservation H . We have 
followed the linear momentum density plot along the lines 
described in M and we have found a behavior similar to 
that of topological solitons at the collision time. As a con- 
clusion, despite that in our case there is no topological 
invariant associated to the qballs, the dynamics related 



to the right angle scattering behavior is the same as that 
for their topological counterparts. We should therefore 
expect such a dynamical behavior to occur generically. 
We shall study in this section the validity of the above 
remarks. 





FIG. 10. Head-on collision of two qballs with the same 
charge. Contour plots of the charge density are given. Con- 
tour levels: 0.4, 0.7, 1, 1.3, 1.6, 1.9, 2.2. Three snapshots are 
shown: initial ansatz (upper entry, time i = 0), at collision 
time (middle entry, t = 41.9) and well after collision (lower 
entry, t = 55.9). Parameter values: B = 4/9, ui = 0.75, initial 
velocity of each qball vo = 0.2. 
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FIG. 11. Head-on collision of two qballs with the same 
charge. Contour plots of the charge density are given. Con- 
tour levels: 0.4, 0.8, 1.2, 1.6, 2.4, 3.2. Snapshots at: t = 
(upper entry), t = 19.6 (middle entry) and t = 69.8 (lower 
entry). Initial velocity of each qball vo = 0.4. Rest of param- 
eters as in Fig. hd. 

Following the simulation for later times, the two qballs 
are seen to eventually stop along the y-axis. They then 
attract, collide again and subsequently scatter and re- 
emerge on the X-axis. This scenario goes on and gives an 
oscillating system with successive right angle scattering. 
We shall not pursue further here this oscillating behavior. 

The next simulation we present has been prepared in 
a way similar to the previous one, but now the initial 



FIG. 12. Head-on collision of two qballs with the same 
charge. Contour plots of the charge density are given. Con- 
tour levels: 0.4, 0.8, 1.6, 2.4, 5, 8. Snapshots at: i = (upper 
entry), t = 12.6 (middle entry) and t = 29.3 (lower entry). 
Initial velocity of each qball vq = 0.8. Rest of parameters as 
in Fig. 

velocity of the qballs is set to a higher value vq — 0.4. 
The results are shown in Fig. |ll|. At the initial stages, 
the simulation does not differ from our first simulation. 
At collision time (middle entry) a central qball is formed. 
However, the subsequent evolution is considerably differ- 
ent. In the lower entry of the figure we have two qballs 
which have evolved from a right angle scattering process. 
They are located on the y-axis and are almost static. In 



addition to that, two other qbaUs are continuing their 
route along the x-axis and they are drifting away from 
each other. Their drift velocity is approximately 0.24. 

It would be interesting to have some insight in the 
above described process. The basic remark is that the 
symmetry that has lead to the right angle scattering in 
the first (slow qballs) case, leads here again to the same 
result. On the other hand, the high kinetic energy of 
the initial qballs allows for more complicated scenarios 
like the formation of new qballs. That is, some amount 
of energy that does not follow the right angle scenario, 
reorganizes to form qballs on the x-axis. In the case 
of topological soliton interactions, arguments related to 
topology would usually preclude such scenarios. These 
arguments do not apply in the case of qballs thus we are 
led to a novel result. An interesting observation is that 
the four solitons which are created after the collision seem 
to be qballs with a frequency different than that of the 
initial ones. This can be seen by comparing the last entry 
of Fig. ^with the profiles of qballs of various frequencies, 
obtained by our numerical shooting method. 

We have repeated the above scattering numerical ex- 
periment with qballs of higher initial velocity vq = 0.8. 
The results of the simulation are given in Fig. U2. A 
new type of collision different than in the two previous 
cases arises. The qballs actually collide to form a central 
qball at the origin (middle entry of the figure). However, 
after the collision (lower entry of the figure) two qballs re- 
emerge traveling along the x-axis and drifting away from 
each other. Their drift velocity after collision is approx- 
imately 0.75. No right angle scattering is present in this 
case so the qballs appear to be almost non-interacting. 
This result is dramatically different from what we see in 
the case of the slow qballs of Fig. ^ where pure right 
angle scattering occurs. It resembles however similar re- 
sults discussed in the literature for scattering of topolog- 
ical solitons moving with very high relativistic velocities 

Notice also that Ward's chiral model, which is inte- 
grable in 2-1-1 dimensions, presents in some respects sim- 
ilar dynamical behavior. Specifically, it has solutions 
which correspond to the right angle scattering S and 
others which represent noninteracting solitons pl| al- 
though the velocity of the solitons does not seem to play 
any role there. On the other hand, the type of scattering 
shown in Fig. O seems to be specific to nontopological 
solitons and, to our knowledge, it has not been observed 
in other two-dimensional models. 

The right angle scattering of solitons has also been 
observed in a three-dimensional Yang-Mills-Higgs theory 
|l2| . The problem has been approached within the mod- 
uli space, that is, the space of the multi-monopole solu- 
tions in the Bogomonly limit. It has been stressed that 
the scattering behavior of the monopoles can be under- 
stood when they are considered almost static during the 
interaction. This requires that the velocity of the inter- 
acting monopoles is small. Such is certainly the case with 
the first of our simulations here shown in Fig hfl. 
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FIG. 13. Head-on collision of two qballs with opposite 
charge. Contour plots of the charge density are given. Solid 
lines represent positive values and dashed-dotted lines neg- 
ative values. Contour levels: ±0.2, ±0.4, ±0.8, ±1.2, ±1.6. 
Snapshots at: t = (upper entry), f = 28 (middle entry) and 
t — 61.4 (lower entry). Initial velocity of each qball vo — 0.4. 
Rest of parameters as in Fig. hd 



Additional insight in the dynamics can be gained by 
looking at collisions between qballs of opposite charge. 
Such are two qballs with the same profile and opposite 
field precession frequency (cf. Eq. (p2|)). Notice that the 
arguments of M can not be applied in this case, at least 



not in a straightforward manner. In particular, these ar- 
guments indicate that we may not expect that the two 
initial qballs can re-emerge after colhsion, travehng at 
a direction perpendicular to the initial one. Rather, we 
would have to think of a combination of parts of the 
initial qballs which would form the final solution after 
collision. Such a combination, involving parts of qballs 
with opposite charge, does not seem to form any solution 
of the present theory, therefore we don't expect a right 
angle scattering behavior in the current case. We have 
to resort to numerical simulations in order to determine 
the actual behavior of the system. 

We now use an initial ansatz appropriately chosen to 
account for opposite-charge qball collision. This ansatz is 
shown in the upper entry of Fig. 03 through contour plots 
of the charge density. The positive-charge qball is located 
on the left of the figure and the negative-charge one on 
the right. We apply to each of the initial qballs a Lorenz 
boost with a moderate velocity vq = 0.4. The middle 
entry of the figure shows that, quite interestingly, at the 
time of collision a mixed state is formed between the two 
qballs. After that, the two qballs separate again, they ap- 
pear to pass through each other and drift apart. The final 
qballs are not identical to the initial ones. The charge Q 
of each of them is approximately half of that of the ini- 
tial state and it corresponds to a frequency lu ^ ±0.6. 
They have also decelerated and their mean velocity is 
approximately equal to 0.25. 

In Sec. II, we have found an interaction between the 
pair of oppositely charged qballs which introduces an os- 
cillation frequency to the system equal to twice the qball 
internal frequency. This interaction is present here and 
it manifests itself in oscillations of the magnitude of the 
field values and also in oscillations of the position of the 
qball centers. The oscillations appear when the qballs 
are coming close to collide. They also persist after the 
collision even when the qballs are well separated {i.e. at 
the instant of the last entry of Fig. |l^) and they don't 
show any tendency to fade out. 

Computer plots of the distribution of the energy den- 
sity of the system shows that some small amount of en- 
ergy is dissipated at right angles after collision. The un- 
derlying mechanism for this phenomenon is presumably 
similar to that leading to right angle scattering in the 
first of our simulations in this series. Similar behavior 
has been observed for two colliding topological solitons 
with opposite topological charge |1^. In this case, the 
solitons annihilate and the energy is dissipated at right 
angles. 

The scattering behavior described by the simulations of 
this section can be viewed as a generalization of the cor- 
responding behavior of topological solitons. We observe 
the two extreme cases which have also been observed in 
the topological case (pure right angle scattering at low 
velocities and pure forward scattering at very high veloc- 
ities) but we also observe an intermediate case of a qball 
split to a forward and a right angle scattering component 
at intermediate velocities. The forward scattering com- 



ponent gets amplified at high collision velocities. This 
amplification occurs at the expense of the right angle 
scattering component. 
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FIG. 14. The fractional charge in forward scattering of 
qballs as a function of the rest frame qball velocity. The for- 
ward component is clearly amplified at large velocities. Points 
in the figure are plotted every 0.25 velocity units. Then the 
points are connected by straight lines. 



The rich behavior of qballs in the scattering simula- 
tions calls for an overview of the different phases. 



• Low velocity scattering of qballs with the same 
charge leads to initial right angle scattering and 
a bound system that performs breather type oscil- 
lations. 



• Intermediate velocity scattering (0.3 < v < 0.7) 
leads to a combination of forward and right an- 
gle scattering and the forward scattered qballs drift 
away from each other and escape to infinity. 



• High velocity scattering {v > 0.7) leads to pure for- 
ward scattering. The forward scattered qballs drift 
away from each other and escape to infinity. 



To clarify the transition from pure right angle scat- 
tering to pure forward scattering in the case of identi- 
cally charged non-topological solitons we plot (Fig. |lj) 
the charge of each of the forward scattered qballs as a 
fraction of the initial qball charge, versus the initial ve- 
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locity ^. Since the total charge is conserved the rest of 
the charge is scattered at right angles. A small amount 
is dissipated throughout our lattice. 

The charge of the outgoing qballs is computed as fol- 
lows: We let the qballs travel 30 space units on the x-axis, 
after the collision and we integrate the charge density 
within a disc of radius 15 space units around each qball 
center. This result is considered to be the charge of the 
final qballs. 

In the case of scattering of oppositely charged qballs 
we have also observed an increase of the forward scat- 
tered component charge for high velocities (upper curve 
of Fig. 14). This increase however does not occur due to 
reduced right angle scattering (this does not happen in 
this case) but due to reduced annihilation between the 
opposite qball charges which is a result of the reduced 
time of overlap between the fast moving qball profiles. 



V. CONCLUSION - OUTLOOK 

We have studied the stability and the dynamics of 
qballs in the context of a simple toy model. We have 
found that the qball instability sector is a very small 
sector in parameter space. The validity of a virial the- 
orem in the stability sector has also been verified and 
we have shown that systems of qball pairs tend to per- 
form breather-type oscillations for long time periods com- 
pared to the period of the internal field rotation. Finally 
we have studied numerically the scattering of qball pairs 
and found that qballs in two space dimensions with the 
same charge tend to scatter at right angles and sub- 
sequently perform oscillations with repeated collisions 
through right angle scattering. At relativistic collision 
velocities a significant forward scattering component has 
also been found which gets amplified as the collision ve- 
locity increases. This effect is consistent with numeri- 
cal experiments involving topological solitons where the 
charge is discretized due to topology. In that case com- 
plete forward scattering is found for v > Vcrit — 0.9 
[Q. In the nontopological case the smooth interpola- 
tion between the above two regimes is allowed because 
the charge is not discretized. 

These results are interesting not only in the context 
of the general study of the solitonic dynamical proper- 
ties but also in the context of realistic physical systems. 
For example in a cosmological setup where dark matter 



comes from susy qballs, breather-type large qball sys- 
tems could lead to detectable signatures in the gravi- 
tational wave spectrum either by direct emission or by 
exciting vibration modes of the neutron stars where they 
can be trapped ||l^. The detectable characteristics of 
such breather type cosmological systems are currently 
under investigation. Another interesting implication is 
related to the statistical mechanics of qball systems [0 . 
Our results imply that the number of qballs is not con- 
served in such systems and in fact it may increase rapidly 
during multiple high velocity collision processes. This is 
in contrast with the assumptions made in recent work |5[| 
which assumes that merging is the only possible outcome 
of qball collisions. 

A natural extension of our work is the numerical study 
of the formation of qballs during a (cosmological) phase 
transition. Such studies based on numerical simulation 
s |lj,|l^ of a system though a temperature quench have 
been performed extensively in the context of topological 
solitons |l^,|l^ but not in the context of nontopologi- 
cal configurations like the ones considered in the present 
study. Mechanisms of thermal creation and annihilation 
of nontopological solitons in the aftermath of a phase 
transition have been previously studied. For small Q 
balls it is assummcd their fusion into larger ones to be 
dominant for both high and low thermal reaction rates 
|l7| Clearly our results modify the above picture in an in- 
teresting way. The numerical simulations performed here 
can be extended in a straightforward way to apply to the 
case of qball formation during a quench. The realization 
of this extension is currently in progress. 



VI. ACKNOWLEDGEMENTS 



We thank N. Papanicolaou for useful conversations 



[1] T.D. Lee and Y. Pang, Phys. Kept. 221 (1992) 251 
[2] S. Coleman, Nucl. Phys. B262 (1985) 263 
[3] A. Kusenko, Phys. Lett. B404 (1997) 285, 



^In order to have clearly separated qballs in the charge com- 
putation. Fig. nil shows results only for v > 0.3 and v > 0.35 
for the opposite and same-charge case respectively. Then 
the two qballs clearly separate from each other after colli- 
sion and drift away. For smaller velocities, the qballs collide, 
go through each other, but do not subsequently drift away. 
Instead the system remains bounded and performs breather- 
type oscillations. 



11 



[4] A. Kusenko, Nucl. Phys. Proc. Suppl. 62A/C (1998) 
248, G. Dvali, A. Kusenko and M. Shaposhnikov, 
Phys.Lett. B417 (1998) 99; K. Enqvist and J. Mc- 
Donald, Phys. Lett. B425 (1998) 309; ibid Nucl. 
Phys. B538, (1999)3221, A. Kusenko, M.Shaposhnikov, 
P.G. Tinyakov and I.G.Tkachev, Phys. Lett. B423, 
A. Kusenko, V. Kuzmin, M. Shaposhnikov and P.G. 
Tinyakov, Phys. Rev. Lett. 80 (1998) 3185, M. Axenides, 
E.G. Floratos and A. Kehagias, Phys. Lett B444 

[5] A. Kusenko and M. Shaposhnikov, Phys. Lett. B418, 
(1998) 104. 

[6] G.H. Derrick, J. Math Phys. 5, 1252 (1964) 

[7] P.J. Ruback, Nucl. Phys. B 296 (1988) 669 

[8] R.S. Ward, Phys. Lett A 208 (1995) 203; 

T. loannidou and W.J. Zakrzewski, J. of Math. Phys. 39 
(1998) 2693 

[9] S. Komineas, "Scattering of solitons in two dimensions" 



cond- mat/9904093. 



[10] E.P.S. Shellard, Nucl. Phys. B 283 (1987) 624 
[11] R.S. Ward, J. Math. Phys. 29 (1988) 386 
[12] N.S. Manton, Phys. Lett B 110 (1982) 54; 

N.J. Hitchin, N.S. Manton and M.K. Murray, Nonlinear- 

ity 8 (1995) 661 



[13] J. Madsen Phys. Lett. B435 (1998) 125, tiep 



ph/980643S 



[14] J. Ye and R.H.Brandenberger, Mod. Phys. Lett. A5 
(1990) 157; Nucl. Phys. B346 (1990) 149 



[15] W.H. Zurek Phys. Rept. 276 (1996) 177, cond- 



mat/9607135 



[16] P. La guna and W.H . Zurek, Phys. Rev. Lett. 78 (1997) 



2519, gr-qc/9607041 



[17] J.A. Frieman, G.B. Gelmini, M. Gleiser and E.W. Kolb, 
Phys. Rev. Letts 60, 2101 (1988); K. Griest and E.W. 
Kolb, Phys. Rev. D40, 3231 (1989); J.A. Frieman, A.V. 
Olinto, M. Gleiser and C. Alcock, Phys. Rev. D40, 3241 
(1989). 



12 



